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Abstract 



A systematic study has been made on the metal-insulator (MI) transition of the doubly 

l o j ; 

degenerate Hubbard model (DHM) in the paramagnetic ground state, by using the slave- 

> 

boson mean-field theory which is equivalent to the Gutzwiller approximation (GA). For 

t^j- ! the case of infinite electron-electron interactions, we obtain the analytic solution, which 
O 

; becomes exact in the limit of infinite spatial dimension. On the contrary, the finite- 



interaction case is investigated by numerical methods using the simple-cubic model with 



the nearest-neighbor hopping. The mass-enhancement factor, Z, is shown to increase 

o 

divergently as one approaches the integer fillings (N = 1, 2, 3), at which the MI transition 

>: 

^ . takes place, N being the total number of electrons. The calculated N dependence of 
Z is compared with the observed specific-heat coefficient, 7, of Sri^La^TiOs which is 
reported to significantly increase as x approaches unity. 
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§1. Introduction 

A study on the metal-insulator (MI) transition has been one of the most fascinat- 
ing subjects in solid state physics. ^ In recent years much attention has been paid to 
transition-metal compounds of A 1 _ :r B a; M03 where A = Sr, Ca, B = La, Y, and M = Ti, 
V, Cr. 2 ' 3 ) These systems exhibit various phases such as a paramagnetic metal (PM), para- 
magnetic insulator (PI) and antiferromagnetic insulator (AFI) when the temperature (T), 
pressure (P), magnetic field (H) or chemical substitution (x) is changed. Theoretically 
the MI transition has been so far investigated by using the single-band Hubbard model 
(SHM), 4 ~ 6 ) with which observed T — P and T — x phase diagrams have been qualitatively 
understood. 7 ~ 12 ) We should, however, remind the fact that constituent transition metals of 
these systems inevitably have the orbital degeneracy. Although the important role played 
by the orbital degeneracy has been recognized, the additional degree of the orbital degen- 
eracy has prevented us from a systematic, theoretical analysis using the degenerate-band 
Hubbard model. 

In the last few years, several theoretical studies have been made of the MI transi- 
tion in the Hubbard model with orbital degeneracy, by employing various methods such 
as the Gutzwiller approximation (GA), 13-16 ) Monte-Carlo method 17-19 ) and slave-boson 
theory. 20 ~ 22 ) By using the GA, Lu 13 ) obtained the critical interaction for the MI transition, 
U c , given by U c ~ (D + 1) at the half filling, while Gunnarsson et al. 17 ' showed U c ~ y/D, 
using a diffusion Monte Carlo method, D denoting the orbital degeneracy. Biinemann 
and Weber 14 ) discussed the first-order MI transition in the half-filled doubly degenerate 
Hubbard model (DHM) when the exchange interaction is included within the GA. The 
MI transition of the infinite-dimensional DHM in the paramagnetic state is investigated 
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by Kotliar and Kajueter, 18 -* and by Rozenberg 19 -* using the dynamical mean-field theory. 

In a previous paper (referred to as I), 20 - 1 the present author proposed the slave- 
boson functional-integral theory for the degenerate-band Hubbard model, by adopting 
the method originally employed for the Anderson lattice model. 23 ) It is a generaliza- 
tion of the slave-boson method of the SHM 24 ) to the degenerate-band model, and its 
mean-field approximation is equivalent to the GA, 13-15 ' 25 ) which becomes exact for the 
Gutzwiller wavefunctions in the limit of the infinite dimension. 16,26 ) An alternative slave- 
boson method was independently proposed by Fresard and Kotliar. 21 ) In our second paper 
(referred to as II) 22 ) the antiferromagnetic state of the half-filled DHM was studied by 
using our slave-boson theory. 20 ) 

Although these studies 13-22 ) are useful in understanding the existing experimental 
data, it is more desirable to make a systematic study of the MI transition of the degenerate- 
band Hubbard model. We have made such a study by using the slave-boson mean-field 
theory developed in I, which is the purpose of the present paper. In the next §2, we briefly 
discuss the mean-field approximation to our slave-boson theory 20 ) applied to the DHM. 
One of the advantages of our slave-boson mean-field theory over sophisticated methods 
such as the dynamical mean-field method, 18 ' 19 ) is that we can obtain analytic solutions in 
some limiting cases. Indeed, we present, in §3, the analytic solution when the electron- 
electron interactions are infinite. For finite interactions, we perform, in §4, numerical 
calculations using the simple-cubic model. The final section (§5) is devoted to conclusion 
and discussion on the experimental data of Srx.^La^TiOs. 3 ) 



2. Formulation 



We adopt the DHM whose Hamiltonian is given by 

H = ^ij C lmcr C jm'<7 + ^ ^mm' C lmCT C i"io-c| m / ./Ci m ' (T ' , (1) 

CT ij mm' i (m,a)^(m' ,a') 

where Q mcr is an annihilation operator of an electron with an orbital index m and spin 
c(— T) I) 011 the lattice site i. The electron hopping is assumed to be allowed only between 
the same sub-band: £™ m ' = tij5 mm > for a simplicity. The on-site interaction, U™^,, is given 
by 

tCm' = U = U iorm = m',a^<7>, (2) 
= \J X = U - 2 J for m ^ m', a ^ a', (3) 
= U 2 = U - 3J for m ^ m', a = a', (4) 

where U and J are Coulomb and exchange interactions, respectively. 

For the paramagnetic state of the DHM, we introduce the boson operators of e*, Pi, 
di, ti and f\ which denote the empty, singly-, doubly-, triply- and fully-occupied states, 
respectively, at a given i site. 20 ' 22 ) As for the doubly occupied states, we take into account 
three kinds of configurations; dio for a pair of electrons on the same orbital with opposite 
spin, dn on the different orbital with opposite spin, and d i2 on the different orbital with 
same spin. These boson operators are under constraints of their completeness and the 
correspondence with fermion operators. 20 ' 22 - 1 

When we apply our slave-boson mean-field theory developed in I to the DHM under 
consideration, we get the ground-state energy given by 

E = E 1 + E 2 , (5) 



where 

Ei = 2qe , 

E 2 = 2(U Q d + U 1 d 1 + U 2 d 2 ) + 2(Uo + U 1 + U 2 )(2t + f). (6) 

Here the occupancies of do, d±, d 2 , t and / are the expectation values of respective boson 
operators and they are determined by the following self-consistent equations: 20,22 -* 

U + e (dq/dd o ) = 0, (7) 

U 1 +e (dq/dd 1 ) = 0, (8) 

U 2 + e (dq/dd 2 ) = 0, (9) 

2(U + Ui + U 2 ) + e (dq/dt) = 0, (10) 

(U + C/i + U 2 ) + e (dq/df) = 0. (11) 



where 



with 
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(Vdo + Vdi + Vd2)(VP + Vt) + Vep + Vt fj\ 2 



N(4-N) ' (12) 

de ^2 e p {e) (13) 

N = 2 f d£ ^ Po{e) (14) 

J cr 

e = l-N + 2(d + d l + d 2 ) + 8t + 3f, (15) 

p = N/4 - (d + d x + d 2 ) -3t-f. (16) 

In eqs. (7)-(16), \x is the Fermi level, po(e) denotes the non-interacting density of states 

per sub-band per spin, and N is the total number of electrons which is 1 and 2 for the 
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quarter- and half-filled bands, respectively. The band-narrowing factor, q, is expressed 
in terms of various occupancies as given by eq. (12). The self-consistent equations given 
by eqs. (7)- (16) are equivalent to those obtained by employing the GA to Gutzwiller's 
wavefunctions. 13-15 ) 

§3. Analytic Solutions for Infinite Interactions 

When the interactions are infinite, we can obtain analytic solutions for the self- 
consistent equations given by eqs. (7)- (16). Solutions for the less-than-half filling (N < 2) 
are classified into the three cases depending on the value of j defined by 



Solutions for the more-than-half filling (N > 2) can be obtained from the result for N < 2 
with the electron-hole symmetry. 

3.1 Case of j = 

When j = 0, three, double occupancies become equivalent: d = d± = rf 2 = d. Taking 
the U — > oo limit of eq. (7), (8) or (9) with 



j = lim (J/U). 



(17) 



t = f = o, 



(18) 



we get 



d = d 



i — 



d 2 = 



(for N < 1), 



(19) 



(1/6) (N-l) 



(for 1< N < 2) 



(20) 



From eqs. (18)-(20), e and p in eqs. (15) and (16) are given by 



e = 1-N 



(for N < 1), 



(21) 



= (for 1 < N < 2). (22) 

p = N/4 (for N < 1), (23) 

= (1/4) (2 — N) (forl<iV<2). (24) 

Substituting eqs. (18)-(24) to eq. (12), we obtain the band-narrowing factor: 

q = 4(1- N)/(4- N) (forJV<l), (25) 

= 6 (N -1)(2- N)/N(4- N) (for 1 < AT < 2). (26) 

Near N — 1 and 2, q is given by 

g = (4/3) 5 (for AT ~ 1), (27) 

= 25 (for ~ 1), (28) 

= (3/2)5 (forA^~2), (29) 

where the electron or hole doping concentration, 5, is given by 5 = \ 1 — A^ | in eqs. (27) 
and (28) and 5 = | 2 - N | in eq. (29). 

3.2 Case of j > 
Similarly we get 

d o = d 1 = t = f = 0, (30) 

d 2 = (for N < 1), (31) 

= (1/2) (N - 1) (for 1< A^ < 2). (32) 

The Af-dependence of e and p is given by eqs. (21)-(24). The band-narrowing factor 
becomes 

q = 4 (1 - N)/(4 - N) (foriV<l), (33) 



= 2 (N -1)(2-N)/N(4-N) (for 1 < iV < 2). (34) 

It is easy to see that q close to N = 1 and 2 is given by 

q = (4/3) 5 (for iV ~ 1), (35) 
= (2/3) 5 (for N ~ 1), (36) 
= (1/2)5 (for TV ~ 2), (37) 

where 5 = \ 1 — N \ in eqs. (35) and (36) and 5 = | 2 — iV | in eq. (37). 

3.3 Case of j < 
We get 27 ) 

= d 2 = t = f = 0, (38) 

d = (for TV < 1), (39) 

= (1/2) (N - 1) (for 1 < iV < 2). (40) 

The iV-dependences of e, p and q are same as in the j > case. 

It is worth to remark that in the SHM, the solution for U — > oo with iV < 1 is given 

by 

e = 1 - N, (41) 

J? = N/2, (42) 

with d = 0, leading to 

q = 2 (1 — N)/(2 — N) (43) 
-25 (for iV ~ 1), (44) 
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where 8 —\ 1 — N\ . We should note that these solutions of the GA for infinite interactions 
discussed above become exact in the limit of the infinite spatial dimension. 16,26 ) Some 
figures depicting the N dependence of the occupancies and the band-narrowing factor for 
the DHM will be shown shortly (dashed curves in Figs. 3-5). 



§4. Numerical Calculations for Finite Interactions 

When the interactions are finite, we cannot obtain the analytic solution for the self- 
consistent equations given by eqs. (7)-(16), which have to be solved by numerical methods 
for given parameters of N, U and J and non- interacting density of states, po(e). We adopt 
the simple-cubic model band, expressed by an approximate, analytic form: 28 ) 



p (e) = A 


9 - uj 2 ' 


1/2 


C 


1 - uj 2 


1/2 




for 


UJ |< 1, 




= A 


9 - uj 2 


1/2 


B 


)-(M- 


-2) 2 ~ 


1/2 


for 1 <| uj |< 3, 




= 




for 


UJ 


|>3, 










(45) 



with uj = e/2t, A/2t = 0.101081, B/2t = 0.128067 and C/2t = 0.02, t being the nearest- 
neighbor hopping. The energy and the interactions are hereafter measured in units of a 
half of the total band width (W/2 = 6t = 1). 



4-1 Results as a function of U 

Firstly we discuss the U dependence of the band-narrowing factor, q, at the integer 
fillings of N — 1 and 2. 29 ) Figure 1(a) shows the U dependence of q for various values of 
J/U of the quarter-filled case (N — 1). When the U is increased, the q factor decreases 
and vanishes above U c which stands for the critical interaction for the MI transition. We 
get q ~ (U c — U) at U < U c , which is the same as in the SHM. 30 ) In all the cases shown 
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in Fig. 1(a), the MI transition is of the second order. The calculated Uc is plotted by the 
solid curve in Fig. 2 as a function of J/U, which shows an increase in U c as increasing 
J/U. On the contrary, the behavior of the MI transition in the half-filled band (N = 2) is 
quite different from that in the N — 1 case, as is shown in Fig. 1(b). The MI transition 
becomes the first- order one for finite J/U, 14 ' 20 ^ and U c decreases as increasing the absolute 
value of J/U. Then the U c curve has a cusp at J/U = 0, as shown by the dashed curve 
in Fig. 2. 

4-2 Results as a function of N 

Figures 3(a) shows the N dependence of the occupancies for various U values with 
J/U = 0. At iV = only the empty occupancy e is finite (e = 1). As increasing N 
from zero, e monotonously decreases while p has a peak at N = 1, where the double 
occupancies, do = di = d 2 = d, begin to increase. Since t and / are one or two-order 
magnitude smaller than the other occupancies, they are not shown in the Figure. Dashed 
curves express analytic solutions of e, p and d for U = oo given by eqs. (18)-(24). We 
note that as the U value is increased, the N dependence of the occupancies approaches 
that for U = oo and that the occupancies for U = 4 are almost the same as those for 
U = oo. 

Similarly, the N dependence of the occupancies for various U values with J/U = 0.1 is 
shown in Fig. 3(b). General behavior for J/U = 0.1 is very similar to that for J/U = 0, 
except that only d 2 has an appreciable magnitude among the three, double occupancies 
because U > U\ > U 2 for J > 0. 

Figures 4(a) and 4(b) show q as a function of iV for various U with J/U = and 
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J/U = 0.1, respectively. The band-narrowing factor is unity at N = 0, and it has dips 
or zeros at N — 1 and 2. Our result of the N dependence of q for J/U = is similar to 
those reported in refs. 19 and 21. The mass-enhancement factor, Z, which is given as the 
inverse of q, for J/U = and J/U = 0.1 is shown in Figs. 5(a) and 5(b), respectively, 
as a function of N. For intermediate interactions (U ~ W), Z has peaks at the integer 
band fillings of N — 1 and 2, where Z shows the divergence for the strong interaction 
(U 3> W). As increasing the U value, Z increases as expected. Dashed curves in Figs. 
5(a) and 5(b) express the results for U = oo given by eqs. (25)-(26) and (33)-(34). We 
note that Z at N < 1 is not symmetric with that at TV > 1 and it is not the same as that 
at N < 2, in particular for J/U = 0.1. This is more clearly seen in the case of U — oo 
as shown by dashed curves in Figs. 5(a) and 5(b) and as given by eqs. (27)-(29) and 
(35)-(37), which show the <5 _1 divergence as approaching the integer fillings. 

§5. Conclusion and Discussion 

By using the slave-boson mean-field theory 20 ) which is equivalent to the GA, 13 ~ 16,25 ) 
we have systematically studied the DHM, which shows the MI transition to take place at 
the integer fillings of N — 1,2 and 3 in the paramagnetic state. When the magnetically 
ordered state is taken into account, this MI transition is modified, particularly in the 
half-filling case, where the antiferromagnetic state can be realized for a relatively small 
interaction. Indeed, for the simple-cubic model with the nearest-neighbor hopping, the 
antiferromagnetic insulator (AFI) is stabilized for an infinitesimally small interaction, and 
then the MI transition does not take place except at the vanishing interaction point in 
the DHM, 22 ) just as in the SHM. 24,31 ) The schematic, phase diagram of the MI transition 
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of the DHM is shown in Fig. 6. When approaching the integer filling from below or above 
with the moderate interaction (U ~ W), the mass-enhancement factor, Z, is considerably 
increased at the PM-PI transition. On the contrary, when the PM-AFI (PM-FI) transition 
occurs with the strong interaction and/or the favorable band filling (as the half filling), the 
mass-enhancement is not realized in the AFI (FI) because of the presence of the energy 
gap in the density of states, 32 ) although the enhancement of Z is expected as approaching 
this transition from the metallic side. Our calculations applying the slave-boson mean- 
field theory to the DHM, 20,22 ) have well accounted for the essential features of the MI 
transition of a system with orbital degeneracy: 

It has been reported that the linear coefficient of the low-temperature specific heat, 
7, and the susceptibility of Sr 1 _ a ,La. E Ti03 significantly increase as x approaches unity. 3 ) 
Filled circles in Figs. 5(a) and 5(b) express the observed mass-enhancement ratio of 
7/7o (— Z). We note that the calculated N dependence of Z well reproduces the ob- 
served x dependence of 7/70. It was previously tried to explain the experimental data 
with the use of the SUM, 12 ) in which the mass-enhancement factor, Z, of its paramagnetic 
state divergently increases near the half filling (eqs. (43) and (44)). We should, however, 
remind the fact that this enhancement does not occur if the antiferromagnetic state is 
included, as was discussed above (Fig. 6). For a better understanding of the MI transition 
in Sr 1 _ a; La a ;Ti03, we have to adopt the triply degenerate Hubbard model in a strict sense, 
because d electrons in this system belong to t 2g bands. 3 ) In order to discuss the high-energy 
phenomena of transition-metal compounds 2,3 ) such as photoemission and the optical con- 
ductivity, it is necessary to take into account " incoherent-quasiparticle" states arising 
from boson fluctuations around the saddle-point solution: 33 ) our slave-boson mean-field 
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theory includes only "coherent-quasiparticle" states, which are relevant to the low-energy 
phenomena like the mass enhancement. Going beyond the mean- field approximation is 
left as one of our future subjects in the slave-boson functional integral theory for the 
degenerate-band Hubbard model. 
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Figure Captions 

Fig. 1 The U dependence of the band-narrowing factor, q, for various J/U with (a) 
N = 1 and (b) N — 2. 

Fig. 2 The critical interaction for the MI transition, U c , in the paramagnetic state as 
a function of J/U for N = 1 and 2, the solid (dashed) curve denoting the second (first) 
order transition. 

Fig. 3 The N dependence of the occupancies for U — 2 (thin solid curve), 4 (bold solid 
curve) and oo (dashed curve) with (a) J/U = and (b) J/U = 0.1. 

Fig. 4 The N dependence of the band narrowing factor, q, for various U with (a) J/U = 
and (b) J/U = 0.1, the dashed curve denoting the result for U — oo. 

Fig. 5 The iV dependence of the mass-enhancement factor, Z(= l/q), for various U with 
(a) J/U = and (b) J/U = 0.1, the dashed curve denoting the result for U — oo. Circles 
express the experimental data of the enhancement ratio of 7/70 in Sr 1 „ a ,La a; Ti03 [ref. 3], 
70 being taken to be the 7 value at N = 0. 

Fig. 6 The electron number vs. interaction phase diagram of DHM (schematic), showing 
the paramagnetic metal (PM), paramagnetic insulator (PI), antiferromagnetic insulator 
(AFI) and ferromagnetic insulator (FI). The relative stability between AFI and FI at the 
quarter filling depends on the electronic and lattice structures and on the interaction of 
a given system. 
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